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ABSTRACT 

We investigate the multiplicity properties of 408 B-type stars observed in the 30 Doradus region of the Large Magellanic Cloud with 
multi-epoch spectroscopy from the VLT-FLAMES Tarantula Survey (VETS). We use a cross-correlation method to estimate relative 
radial velocities from the helium and metal absorption lines for each of our targets. Objects with significant radial-velocity variations 
(and with an amplitude larger than 16 kms“^) are classified as spectroscopic binaries. We find an observed spectroscopic binary 
fraction (defined by periods of < 10^'® d and mass ratios > 0.1) for the B-type stars, fsiobs) — 0.25 ± 0.02, which appears constant 
across the field of view, except for the two older clusters (Hodge 301 and SL 639). These two clusters have significantly lower binary 
fractions of 0.08 ±0.08 and 0.10 ±0.09, respectively. Using synthetic populations and a model of our observed epochs and their 
potential biases, we constrain the intrinsic multiplicity properties of the dwarf and giant (i.e. relatively unevolved) B-type stars in 
30 Dor. We obtain a present-day binary fraction /B(true) = 0.58 ± 0.11, with a flat period distribution. Within the uncertainties, the 
multiplicity properties of the B-type stars agree with those for the O stars in 30 Dor from the VETS. 
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1. Introduction 


The presence of a close companion can significantly influence 


the evolution of massive stars (e. g. Podsiadlowski et aHIl992t 
iLanger et alJl200^ lEldridge et al.l 201 lb . but the prevalence of 


mass transfer in such systems has only recently been realised. 


ferred for massive stars in the Galaxy (e.g. [Mason et al.l 

2009 

o 

o 

o 

o 

o 

c 

Kiminki & Kobulnickv 

2012 

Chini et al.l2012;lSota et al.l2014 

) is that the majority of massive 


stars may be born as part of, and reside with in, multiple systems 
dKobulnickv & Frven2007HSana et alJ2012h . Additionally, stars 
that curr ently appear to be sin gle may have evolved from binary 
systems dde Mink et al.ll2()T4l) . with important consequences o n 
their inferred evolutionary histories dde Mink et al.ll20fTll2013h . 

The multiplicity fraction of B-t ype stars has been investi¬ 
gated in some Galactic clusters (e.g. lRaboudll 19961) but, to date, 
most studies of the multiplicity of high-m ass stars have focused 
on the more massive O-type objects (e.g. ISana & Evan^l201 iL 
and references therein). Substantial spectroscopic samples of B- 
type stars in the Galaxy and Magellanic Cl ouds are available 
from large spectroscopic p r ogram mes (e.g. lEvans et^ 1200.51 
120061: iMartavan et alJ[200^ l2007l) . but they have lacked suffi- 
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dent time sampling to investigate the role of binarity in the ob¬ 
served populations. 

One of the primary motivations for the VL T-FLAM ES 
Tarantula Survey (VFTS. lEvans et al.]l201 iL here after [Paper ih is 
to investigate the multiplicity properties of the massive stars in 
the 30 Doradus region of the Large Magellanic Cloud (EMC). 
To this end, the VETS obtained multi-epoch spectroscopy of 
^800 O- and early B-type stars. The observed spectroscopic 
binary fraction, /o(obs), of the 360 O-type stars in t he VETS 
was found to be /o(obs) = 35 ± 3% dSana et al.ll2013L hereafter 
iPaper VIIll) : once the observational biases were taken into ac¬ 
count, this gave an intrinsic binary fraction, /o(true) = 51 ± 4%. 


Here we investigate the multiplicity properties of the B-type 
stars observed by the VETS. One of our objectives is an esti¬ 
mate of their present-day observed spectroscopic binary fraction, 
/B(obs), where / defines the fraction of targets where Doppler 
shifts are detected (rather than the total fraction of all stars), 
probing separations with orbit al periods of up to ap proximately 
10 years (e.g. see Fig. 1 from Sana_& Eyand HOllli . For consis- 
tency with the analysis bv ISana et al.l d2013 and in iPaper VIIlL 
we adopt the same definition of a spectroscopic binary as a sys¬ 
tem with a period of <10^ ® d and a mass ratio of >0.1. The 
present-day binary fraction can then be used to constrain the 
minimum binary fraction of B-type stars at their birth (with ini- 
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tial masses of > 8 Mq). Ultim ately, in com bination with the re¬ 
sults for the O-type stars tfrom ll^er Vlllh . our aim is to deter¬ 
mine the fraction of massive stars that might be affected by bi¬ 
nary interaction during their lives, which will influence their evo¬ 
lution and the properties of their final explosions . Indeed, assum - 
ing the form of the stellar mass-function from iKrounal (l2001ll . 
around two thirds of the progenitors of core-collapse supernovae 
were probably early B-type stars when on the main sequence. 

After a brief summary of the observations in Section |2] we 
describe the methods used to identify spectroscopic binaries in 
Section [3] We present our results in Section |4] and investigate 
the intrinsic properties of the binary population in Section |5] 
Concluding remarks are given in Section|6l 

2. Observations 

The VFTS spectra were obtained with t he Fibre Large Array 
Multi-Element Spectrograph (FLAMES, iPasauini et al.l 120021) 
on the Very Large Telescope (VET). Classifications for the 438 
B-type stars observed in the VETS were given by lEvans et aP 
( l2015h . These are located in the main clusters in the 30 Dor re¬ 
gion (i.e. NGC 2070, NGC 2060, Hod p301, SL6 39) and the 
local field population (see Eig. 4 from lEvans et al.h . All of the 
spectra were obtained using the fibre-fed Medusa-Giraffe mode 
of ELAMES, so the sample does not include stars in R136, the 
young massive cluster at the core of 30 Dor, which is too densely 
populated for effective use of the Medusa fibres. 

This paper presents a radial-velocity (RV) analysis of the 
multiple observations of the B-type stars with the LR02 setting 
of the Giraffe spectrograpfQ. The observations spanned 10-12 
months for most targets, with an extended baseline of 22 months 
for 31 targets, due to a reobservation for operational reasons. 
Two or three consecutive exposures were obtained on any given 
night, with the longer-term sampling designed to optimise detec¬ 
tion of binaries with periods up to ^200 d, given scheduling con¬ 
straints. (The detection probability de creases sign ificantly for 
longer-period systems; see e.g. Eig. 8 of lPaner VIIll) . Eurther de¬ 
tails regarding target s election , data reduction, and the observa¬ 
tion dates are given in lPaner j 

Absolut e RVs for the B-type sample were provided by 
lEvans et aH (12015h . which are either the mean of estimates from 
each observation or, where binary motion was suspected, the es¬ 
timate from a particular observation. Eor 30 stars the spectra 
were not of sufficient quality to allow reliable estimates of the 
absolute RVs, and these stars have been excluded from our sam¬ 
ple. As these span a range of spectral types and morphologies 
(including examples of Be-type stars), they should not unduly in¬ 
fluence the results. A further five targe ts are the double-lined bi¬ 
naries identified ^|E^ns|et^ (1203^ namely: VETS 112, 199, 
240, 637, and 69® It was challenging to obtain multi-epoch RV 
estimates to the same standard as those used for the analysis in 
Section [2 but their binary natures were easily revealed by vari¬ 
ations in their spectra (indicative of binary motion, rather than 
simply composite spectra). 

* The LR02 observations provided coverage of 3960 to 4560 A at 
a spectral resolving power of ~6 500. Observations were also obtained 
using the LR03 and HR15N settings, but with considerably limited time 
sampling (cf. the LR02 data) so these were not considered further here. 

^ A detailed discussion of VFTS 698 was given by lDunstallet_al, 
(l2012h . The SB2 system VFTS 652 was omitted from the lEvans et al, 
sample, and henc e from this paper, be cause of the O-type classification 
of the secondary dWa lbom et al.ll2014l) . The ten BO-type stars classi fied 
bv lWalborn et all ( 120141) were also discussed by lEvans et all )l2015h . so 
are included in this paper. 


Therefore, in total we investigated the multiplicity of 408 
stars, comprised of 361 dwarfs and giants (i.e. relatively un¬ 
evolved objects) and 47 supergiant s (defined as having a loga¬ 
rithmic gravity, logg<3.3dex, see IMcEvov et aU 120151) : here¬ 
after these are referred to as the ‘unevolved’ and ‘supergiant’ 
samples. 

3. Radial velocities 

3.1. Methodology 

We began by examining the 403 single-lined targets for evi¬ 
dence of RV variability. An important aspect of RV studies in 
early-type stars concerns the choice of diagnostic lines. Our pre¬ 
ferred features for B-type stars would be their intrinsically nar¬ 
row metal absorption lines. However, the large variation in the 
intensity of these lines with spectral type means they would 
not provide a homogeneous set of estimates on their own. We 
therefore employed both helium and metal absorption lines (see 
Table [T]( to obtain a range of independent RV estimates, as per¬ 
mitted by the signal-to-noise (S/N) ratio of the data and the spec¬ 
tral type of each target. 

For each of the lines available from Table[T] relative RV esti¬ 
mates (vi) were obtained for each LR02 spectrum (*) of a given 
target from cross-correlation of the data with an adopted ref¬ 
erence spectrum (typically that for the object with the largest 
number of counts in the continuum arou nd 4200 A). This ap¬ 
proach has been widely used (see e.g. lMaroccoet21120151 : 
iManick et al.l 120151: iRobertson & Mahadevanll20l4 and refer- 
en ces therein), wi t h our methodology being similar to that used 
bv iMarocco et al.l (l2015l) . 

The spectral region used in the cross-correlation analysis de¬ 
pends on the width of the spectral feature being considered. In 
turn this depends on both its intrinsic width (e.g. diffuse he¬ 
lium lines will be far broader than non-diffuse helium lines 
or metal lines) and on the stellar projected rotati onal velocity. 
The p rojected rotational velocity estimates from iDufton et al.l 
(12013h for the subset of our sample which have no significant 
radial-velocity shifts cover a wide range. The distribution is bi- 
modal with peaks at Ve sin 20 and 180 km s“^, coupled with 
a long tail to around 400 kms“^. The latter would equate to 
a full-width-half maximum rotational broadening of about 6A 
(depending on the stellar wavelength). Given this wide varia¬ 
tion in the width of spectral features, the choice of spectral re¬ 
gions for each line was made subjectively. The principle criteria 
were to minimise the region selected (to maximise the cross¬ 
correlation signal), whilst ensuring both complete coverage of 
the feature and sufficient continuum to fully characterise the 
cross-correlation function. 

The position of the maximum could be affected by spuri¬ 
ous structure in the cross-correlation function, so our method 
was informed by the pr ojected equatorial ro tation velocities 
(Ve sin i) for each star from lPufton et al.l (l2013h . For supergiants 
and narrow-lined spectra (ue sini < 150kms“^) the cross¬ 
correlation function was smoothed by a 5-point moving average 
(equivalent to ~50kms“^ in velocity space). For stars rotating 
more rapidly, a 55-point moving average was used, to militate 
against the increased influence of random noise in the cross¬ 
correlation function; this method was tested for the narrow-lined 
spectra giving, in general, identical results. 

The procedure is illustrated for two cases in Fig. [T] viz. the 
Si III 4553 A line in the narrow-lined bright supergiant VFTS 003 
(B1 Ia+) and the He 14388A diffuse line in the relatively rapidly 
rotating (ve sini~ 212 kms“^) and fainter main-sequence star. 
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Fig. 1. Examples of the methodology to estimate radial-velocity shifts between the spectra. Left-hand panels: VETS 003 is a rel¬ 
atively bright narrow-lined supergiant. The upper panel shows the Si III 4553A line in the template spectrum (MJD = 54804.148, 
black line) and at MJD = 55108.309 (red). The lower panel shows the original (black line) and smoothed (red) cross-correlation 
function in velocity space; no significant velocity shifts were found for this target. Right-hand panels: VETS 189 is a relatively faint, 
broad-lined, main-sequence star. The upper panel shows the He I 4388A line in the template spectrum (MJD = 54867.160, black) 
and at MJD = 54824.242 (red). The lower panel shows the original and smoothed cross-correlation functions in velocity space; the 
estimated radial-velocity shift between these was ~200kms“^. 


Table 1. Absorption lines used to estimate relative radial veloc¬ 
ities for the B-type spectra. 


Ion 

A (A) 

Type 

On 

4080 

Group 

He I 

4009 

Singlet 

He I 

4026 

Triplet 

He I 

4121 

Triplet 

He I 

4144 

Singlet 

He I 

4169 

Singlet 

On 

4317 

Doublet 

He I 

4388 

Singlet 

On 

4414 

Doublet 

He I 

4438 

Singlet 

He I 

4471 

Triplet 

Sim 

4553 

Singlet 


Notes. The OII A4080 group comprises a number of lines extending 
over the 4069-4095 A region. 


VETS 189 (BO.7: V); as such they illustrate the range in quality 
of the available observational data. 

The relatively narrow absorption lines and high S/N of the 
supergiant spectra meant that all of the metal and helium lines 
listed in Table[T]could be used for the RV analysis; this was also 
possible for ~20% of the (non-supergiant) narrow-lined objects. 
Eor more-rapidly-rotating stars {Ve sin *>150 km s “ ^), only the 
helium lines were generally available. 

The relative RVs for each observation of 370 of our targets 
are presented in Table |2l the remaining 33 targets were affected 


by nebular contamination and are discussed in the next section, 
while a RV analysis was not attempted for the five targets classi¬ 
fied as SB2 systems. We note that non-zero estimates arise from 
cross-correlation of the template spectra with themselves, with a 
systematic median offset of —0.28kms“^ arising from the fits 
to the cross-correlation function. In the context of our analysis 
of relative velocities and the adopted criteria to identify binaries 
(Section l33T l such measurement uncertainties are not significant. 

3.2. Nebular contamination 

There is considerable nebulosity across the 30 Dor region, which 
contaminates the EL AMES spectroscopy to varying degrees. 
This can complicate the RV estimates by the introduction of 
additional structure (centred at ARV ~0kms“^) in the cross¬ 
correlation profiles of the helium lines. Eor the majority of our 
targets the nebular emission was sufficiently weak that any as¬ 
sociated structure in the cross-correlation function was minor, 
and robust RV estimates were possible. However, the contamina¬ 
tion was more significant for 33 of our targets, so we employed 
a modified method where all exposures for a given date were 
merged to improve the S/N, and then analysed as for the other 
spectra. The relative velocities from this approach for these tar¬ 
gets are given in Table [3] 

3.3. Criteria for binarity 

Our strategy for binary detection relies on statistical criteria, 
with the signihcance threshold adopted such that it limited the 
number of false positives. Thus, while some binaries may re¬ 
main undetected, this systematic approach enables modelling of 
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Table 2. Relative radial-velocity (RV) estimates for the available lines in each observation (cf. the adopted template spectrum for 
each star). Values for the hrst three VFTS targets are shown; the full table is available online. 


VFTS 

MJD 





Relative radial velocities [km s ^ ] 







Si III 
A4553 

Oil 

A4080 

Oil 

A4317 

Oil 

A4414 

He I He I He I He I 

A4026 A4471 A4009 A4144 

He I 
A4388 

He I 
A4121 

He I 
A4438 

He I 
A4169 


001 

54815.254 





-25.48 




10.55 




001 

54815.277 





-6.33 




-31.42 




001 

54815.305 





-13.05 




-1.63 




001 

54815.328 





-22.53 




-28.72 




001 

54859.211 





-9.27 




-9.76 




001 

54859.230 





6.93 




-19.24 




001 

54891.062 





-32.84 




-13.82 




001* (S) 

54891.082 





-0.44 




-0.28 




003 

54804.082 

0.76 

2.83 

-0.57 

0.92 

-0.44 

1.05 

-0.05 

-1.98 

2.43 

-0.56 

-1.86 


003 

54804.105 

0.76 

1.33 

0.79 

0.92 

-1.91 

1.05 

1.44 

-0.55 

1.08 

0.88 

0.81 


003 

54804.125 

0.76 

2.83 

0.79 

0.92 

1.04 

1.05 

1.44 

-0.55 

1.08 

0.88 

0.81 


003* (V) 

54804.148 

-0.54 

-0.17 

-0.57 

-0.42 

-0.44 

-0.28 

-0.05 

-0.55 

-0.28 

-0.56 

-0.53 


003 

54804.168 

0.76 

1.33 

0.79 

-0.42 

-0.44 

1.05 

-0.05 

-1.98 

-0.28 

-0.56 

-0.53 


003 

54804.191 

0.76 

-0.17 

0.79 

-0.42 

1.04 

1.05 

1.44 

-0.55 

-0.28 

-0.56 

-0.53 


003 

54836.219 

4.65 

2.83 

2.16 

0.92 

8.40 

10.35 

2.92 

2.30 

3.78 

2.31 

6.14 


003 

54836.238 

4.65 

2.83 

2.16 

0.92 

8.40 

10.35 

4.41 

2.30 

5.14 

2.31 

0.81 


003 

54836.266 

3.35 

2.83 

2.16 

0.92 

8.40 

10.35 

4.41 

2.30 

5.14 

2.31 

0.81 


003 

54836.285 

4.65 

2.83 

2.16 

0.92 

6.93 

11.67 

2.92 

2.30 

3.78 

2.31 

3.47 


003 

54867.086 

13.73 

8.84 

11.72 

7.64 

17.24 

19.64 

13.33 

12.28 

14.62 

10.92 

14.14 


003 

54867.109 

15.02 

10.34 

13.09 

6.29 

17.24 

20.97 

13.33 

12.28 

13.26 

8.05 

7.47 


003 

55108.309 

-0.54 

-0.17 

-0.57 

-0.42 

-0.44 

-1.61 

1.44 

-1.98 

-2.99 

0.88 

-3.19 


005 

54815.254 





-7.80 



6.58 

-1.63 




005 

54815.277 





-0.44 



-0.55 

2.43 




005 

54815.305 





18.71 



-17.67 

3.78 




005 

54815.328 





-7.80 



16.56 

-1.63 




005 

54859.211 





-1.91 



17.99 

24.10 




005 

54859.230 





-10.75 



8.00 

13.26 




005 

54891.062 





-12.22 



-6.26 

-1.63 




005* (S) 

54891.082 





-0.44 



-0.55 

-0.28 




005 

55113.309 





-13.69 



0.87 

21.39 




005 

55113.328 





-0.44 



10.86 

11.91 





Notes. Celestial coordinates and optical photometry for the VFTS stars are given in Table 5 of Paper I. Observations indicated by the * qualifier 
to the first column were the template spectra for the RV estimates. The subsequent letters in parentheses indicate the binary status from this work: 
(B) = binary; (S) = single; (V) = RV variable; see Section [T3] for further details. 


the observational biases. For consistency, we adopted similar cri- 
teria for bin arity to those used in the analysis of the O-type stars 
(iPaper VIlJ) . such that an object is considered a spectroscopic 
binary if at least one pair of RV estimates satishes simultane¬ 
ously 


-!tLjiL>4,0 and \v, 


Vj \ > ARV, 


( 1 ) 


where Vi is the mean relative RV for exposure i with respect 
to the template, and ai is its associated standard deviation. The 
choice of 4.0 for the conhdence threshold for a detection was 
guided by the number of false positives one may statistically 
expect given the sample size, the uncertainties of the RV esti¬ 
mates, and the number of RV pairs. Using Monte Carlo simula¬ 
tions, adopting 4.0 leads to fewer than 0.2 false detections for the 
whole sample, which is sufficient for our purpose s. Furtherm ore, 
this threshold value is identical to that adopted in iPaper VIIlL en¬ 
suring a consistent approach. 

The choice of ARVmin was particularly important for the 
supergiants, for which precise RVs could be determined given 
their high S/N and relatively narrow-lined prohles. Small ranges 
of variations (~5-10kms“^) could indicate a binary compan¬ 
ion, but may also resul t from atmospheric activity or pulsations 
(e.g. lTavlor et al.ll2014l) . 


The percentage of binary candidates identihed for different 
values of ARVmin is illustrated in Fig. For the unevolved 
stars, no clear break is seen. However, the fraction of RV- 
variable objects for the supergiants presents a clear kink around 
ARVniin = 16kms“^ with samples plausibly dominated by bi¬ 
narity and intrinsic variability above and below the kink, re- 
sp ectively. Th is is a slightly lower threshold than that adopted 
in I Paper VIIll for the O-type stars (i.e. ARVmin =20kms“^), 
but is consistent with the lower masses of early B-type stars (cf. 
those for O star^. 

Objects fulhlling both criteria of Eq. [T] (with 
ARVniin = 16kms“^) are considered as spectroscopic bi¬ 
naries. Those fulfilling the first condition of Eq. [T] but with a 
maximum ARV in the range of 5 to 16kms“^ are identihed as 
‘RV variables’, with the variations attributable to either binarity 
or intrinsic variability. The remaining objects are presumed to 
be single stars for the purposes of the following discussionQ 


^ For a given mass ratio and orbital period, the semi-amplitude of the 
RV signal for the primary is oc . 

As discussed in Sectionj^ 30 targets were excluded from our sam¬ 
ple because absolute RV values could not be measured. However, rela¬ 
tive velocity estimates (often only from one feature) were determimed 
for 15 of these targets, with tentative evidence of RV variations found in 
seven stars: VFTS 155, 377, 391 434, 442, 742, and 821 (and flagged as 
such bv lEvans et alJ20f^ . with no significant variations seen in the re- 
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^RV^in (km s ') 

Fig. 2. Upper panel: Inferred binary fraction as a function of 
the adopted threshold velocity (ARVmin) for the unevolved 
(dwarf/giant) B-type stars in our sample (solid line). The dashed 
line shows the corresponding fraction that would have been 
classified as binary candidates but failed the ARVnUn criterion. 
Lower panel: As for the upper panel, but for the supergiants. In 
both plots the vertical dotted line indicates the adopted threshold 
of 16kms“^. 


4. The binary sample 

4.1. Observed binary fraction 

Employing the criteria in Eq. [T] to analyse the RV estimates in 
Tables |2] and [2 combined with the SB2 systems, we found 90 
binaries in the unevolved sample and 11 binaries in the super¬ 
giants. The observed binary fractions are therefore 25 ± 2% and 
23 ± 6%, respectively, where b inomial statistics have been used 
to compute the erTor bars (see lSana et al.ll2()()^ . An additional 
23 unevolved stars and 17 supergiants were detected as RV vari¬ 
ables, comprising 6 ± 1% and 36 ± 7% of their respective sam¬ 
ples. The lower limit for the spectroscopic binary fraction for the 
entire VETS B-type sample {including the supergiants and SB2 
systems) is /B(obs) = 25 ±2%, with RV variables accounting 
for a further 10 ± 1% of the stars. 


maining eight stars: VFTS 366, 407, 408, 462, 573, 653, 689, and 854. 
Given the inhomogeneity of these estimates in terms of the lines and 
observations these stars are not conside red further, but we note these 
results for completeness icf. lEvans et all) . 


Table 3. Relative radial-velocities (RV) for spectra with signif¬ 
icant nebular contamination but which estimates were possible 
from combining the spectra from each night (see Section lO l. 


Values for the first three entries are shown; the full table is avail¬ 
able online. 

VFTS 

MJD 

Relative radial velocities [km s 




He I 
A4388 

He I 
A4144 

He I 
A4026 

He I 
A4471 

004* (S) 

54767 

0.39 

0.44 

0.29 

-1.42 

004 

54827 

-14.28 

-15.10 

-12.47 - 

-71.61 

004 

54828 

-39.97 

-59.7 

-23.65 - 

-14.02 

004 

54860 

-10.62 

-3.45 

17.85 

0.37 

004 

54886 

-60.14 

-40.36 

-12.47 - 

-23.02 

004 

55114 

53.58 

-3.45 

-58.75 - 

-28.42 

010 

54767 

26.80 

5.19 

-14.43 

6.94 

010 

54827 

9.20 

-12.24 

-26.21 - 

-32.39 

010* (S) 

54860 

-0.28 

0.44 

-1.17 

-0.34 

010 

54886 

-39.55 

40.07 

-8.54 


010 

55114 

4.88 

41.65 

-17.38 


054* (S) 

54767 

-0.28 

-0.27 

6.93 

-0.28 

054 

54827 

20.03 

-28.95 

-15.16 - 

-18.88 

054 

54860 

49.82 

24.11 

8.40 

19.64 

054 

54886 

-21.94 

-23.22 

-60.82 - 

-37.47 

054 

55114 

-21.94 

31.28 

-65.24 - 

-33.49 


Notes. Celestial coordinates and optical photometry for the VFTS stars 
are given in Table 5 of Paper I. Observations indicated by the * qualifier 
to the first column were the template spectra for the RV estimates. The 
subsequent letters in parentheses indicate the binary status from this 
work: (B) = binary; (S) = single; (V) = RV variable; see Section lT3] for 
further details. 


4.2. Supergiants 

The effects of macroturbulence in b lue supergiants appear t o 
be linked to line-profile variations (ISimon-Dlaz et al.l 1201 Ol). 
so adoption of ARV^in = 16kms“^ in Section lT3] was par¬ 
ticularly important for these objects in the VFTS sample (for 
which the uncertainties on the RV estimates are relatively small). 
To investigate the nature of RV variations in the supergiants 
further we employed th e same procedures as those used by 
ISimon-Dlaz et al.l (1201 Oh to quantify the potential role of macro¬ 
turbulence in our targets. 

In general, this additional analysis confirmed the classifica¬ 
tions regarding multiplicity using the criteria from Section 13.31 
However, there are four objects (VFTS 420, 423, 541, and 
829) classified as RV variables (with 12 < ARV < 15kms“^) 
for which the results are consistent with minimal macroturbu¬ 
lence, thereby strengthening the case for orbital RV variations. 
Conversely, VFTS 591, classified as a binary using the above 
criteria, has a relatively large contribution from macroturbu¬ 
lence and thus might be a single star. Revising these classifica¬ 
tions would increase the binary fraction for the supergiants from 
23 ±6% to 30 ±7%, i.e. within the statistical uncertainties, in¬ 
dicating that uncertainties due to macroturbulence do not appear 
to strongly influence the result. 

4.3. Timescaies for detected variations 

The upper panel of Figure |3] shows the cumulative distribu¬ 
tion of the maximum RV variation for our objects with es¬ 
timates in Tables |2] and |3] which are classified as binaries. 
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Fig. 3. Cumulative fractions of the maximum RV amplitude de¬ 
tected (upper panel) and the minimum time difference corre¬ 
sponding to significant RV-variations (lower panel) for the de¬ 
tected binaries in the B-type sample. 



Ks mag 



Ks mog 


Fig. 4. Observed binary fraction as a function of magnitude. 


4.5. Brightness variations 

We also investigated the binary fraction as a function of the 
brightness of our targets, employing if^s- band magnitudes from 
the VISTA Magellanic Clouds Survey (ICioni et al.l l201lll . As 
shown in Fig. |4] there are no obvious trends, suggesting that the 
binary fraction is mostly uniform over more than two orders of 
magnitude in brightness. 


Approximately ^20% of the objects in the binary sample have 
ARV > 100km s“^; these are probably short-period binaries. 

The lower panel of Figure |3] shows the minimum time dif¬ 
ference for which significant RV-variations are detected. The 
minimum elapsed time between the LR02 observations of eight 
of the nine fibre configurations was typically between one and 
seven days (except for Field C, with a minimum separation of 
28 days, see Appendix of Paper I); hence we expect that the RV 
detectability is principally determined by the properties of the 
stars rather than the cadence of the observations. Variations are 
seen over timescales of one to ten days for ~30% of the sample 
(although this should not be taken as a direct measure of their 
orbital periods), while a relatively small number of the detected 
systems (~15%) only show variations over several hundreds of 
days. Full characterisation of the individual orbital properties 
of the detected binaries will require comprehensive multi-epoch 
spectroscopy. 


4.4. Spatiai variations 


We searched for variations in the observed binary fraction in 
the different structures in the region, i.e. NGC 2070, NGC 2060, 
H odge 301, SL639 ( adopting the same definitions as those used 
by If vans et al.ll20T5 ). and the field population. The binary frac¬ 
tions and numbers of targets in each subsample are summarised 
in Table |4] The observed fractions are marginally larger in 
NGC 2070 and NGC 2060, but remain within the estimated un¬ 
certainties for those in the field population. For the two older 
clusters, the observed fractions for the unevolved targets appear 
to be slightly lower. 


5. Intrinsic muitipiicity properties 

The observed binary fraction of the 408 B-type stars is 
/B(obs) = 0.25 (i.e. 96 positive detections through RV variations 
and 5 SB2 system). Because of the cadence of the FLAMES 
observations, the S/N of our data, and the orientation of the bi¬ 
nary systems with respect to our line of sight, some binaries will 
have eluded detection. Thus, ignoring statistical uncertainties, 
/B(obs) represents a lower limit on the true binary fraction of 
B-type stars in 30 Dor. 

To estimate the intrinsic binary fraction we modelled our 
criteria for binarity and estimated our observational biases and 
detection sensitivities. In this analysis we ignored the five SB2 
systems discussed in Section |2] as their double-lined nature pre¬ 
vented us from obtaining reliable estimates of the RV variation 
of the primary stars. Therefore, as a result of this approximation, 
the intrinsic binary fraction that we recover may be lower than 
the true value by a few percentage points. However, this remains 
well within the unce rtainty of th e method (see below). 

As discussed in IPaner VIlS the detection probability of a 
given binary system depends on the properties of the data (un¬ 
certainty of RV estimates, time sampling) and on the proper¬ 
ties of the RV signal that one is trying to detect. The latter 
is directly determined by the orbital properties of the system. 
The (in)completeness of the VFTS campaign therefore depends 
on the distribution of the orbital parameters of the parent bi¬ 
nary population (predominantly the orbital periods, mass ratios, 
and eccentricities). Unfortunately, these quantities are unknown. 
Thus, to constrain the intrinsic multiplicity properties of the bi- 
naries in o ur sample, we follow the approach implemented in 
iPaner VIl3 for the analysis of the O-type stars. 
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Table 4. Observed binary fraction for various subpopulations in the 30 Dor region. 


Region 

Non-supergiants 
Nohj /B(obs) 

Supergiants 

Nohj /B(obs) 

Aobj 

All 

/B(obs) 

NGC 2070 

93 

0.27 ± 0.05 

16 

0.31 ±0.12 

109 

0.28 ± 0.04 

NGC2060 

29 

0.34 ±0.05 

3 

0.00 ±0.11 

32 

0.31 ±0.08 

Hodge 301 

12 

0.08 ±0.08 

3 

0.00 ±0.27 

15 

0.07 ± 0.06 

SL639 

10 

0.10 ±0.09 

4 

0.25 ±0.27 

14 

0.14 ±0.09 

Field 

217 

0.24 ± 0.03 

21 

0.24 ±0.22 

238 

0.25 ±0.03 

All 

361 

0.25 ±0.02 

47 

0.23 ±0.06 

408 

0.25 ± 0.02 


Using a Monte Carlo method, we synthesised populations of 
stars with specihed parent orbital distributions. Taking into ac¬ 
count the measurement uncertainties and the time sampling of 
each object, we sought to reproduce the observational quantities 
in three aspects: (i) the observed binary fraction /B(obs), (ii) the 
cumulative distribution of the maximum amplitude, CDF( ARV), 
of significant RV variations (see Fig. [3j upper panel) and (iii) 
the cumulative distribution of the minimum time separation, 
CDF(AHJD), between any pair of RV points that differ signifi¬ 
cantly from one another (see Fig. [3] lower panel). 

The quality of the match between observations and simula¬ 
tions was assessed using a global merit function (S')- This was 
dehned as the product of the Kolmogorov-Smirnov (KS) prob¬ 
abilities between the synthetic CDFs for ARV and AHJD and 
their observed distributions, and of the Binomial probability that 
describes the chance to obtain the same number of detected bina¬ 
ries (A^bin) as in our sample, given the simulated observed binary 
fraction /ejobs)®™”^ and the sample size N: 

S' = Pks(ARV) X Pks(AHJD) x B{Nun, (2) 

The method is described further in lPaner VIIll Stellar masses 
in the simulated populations were randomly drawn from a mass- 
function over a given stellar range. We can reasonably assume 
that the unev olved stars follow a standard mass-function (e.g. 
lKroupall2001l) . but the situation may be different for the super¬ 
giants, many of which would have been born as higher-mass 
stars, which could have experienced considerable evolution prior 
to their current phase (e.g. mass lost via stellar winds, binary in¬ 
teractions). Therefore, we focus on the 357 unevolved B-type 
stars in our sample (which does not include the four SB2 sys¬ 
tems with unevolved stars; cf. the numbers in Table |4]i. 

As in iPaner VlTS we considered orbital periods (P), mass- 
ratios {q = M 2 /M 1 ) and eccentricities (e) over ranges of 
0.15 < log(P/d) <3.5, 0.1<q<1.0 and 0<e<0.9, respec¬ 
tively. We adopted power-laws to describe the distributions of 
the orbital parameters: /p ex (log P/d)’', fqocq'^ and feoce^. 
We varied tt, k and the intrinsic number of binaries to best re¬ 
produce the observed quantities CDF(ARV), CDF(AHJD) and 
/B(obs). Our method was insensitive to the eccentricity distri¬ 
bution and we could not constrain p. As in our analysis of the 
O-stars, we adopted the value ob tained from O-typ e binaries in 
Galactic open clusters, p = —0.5 (ISana et al.ll2012l) . While there 
is no guarantee that this value is appropriate, it preserves the ho¬ 
mogeneity of the analysis of the O- and B-type populations from 
the VFTS. The range of values explored for /B(true)> ^ k are 
given in Table |2 and the criteria to identify binarity in the Monte 
Carlo simulations are those described in Section 

Fig.|5] shows the behaviour of 5' when varying /B(true)> 
TT and K and Fig. |6] shows the simulated distributions for the 
optimum 5' compared with the observations. The best repre¬ 
sentation of the data is obtained with fBrtmpt =0.58 ±0.11, 
TT = 0.0 ± 0.5 and k = -2.8 ± 0.8. 


We estimated the uncertainties on the retrieved multiplicity 
properties by the use of 50 synthetic datasets which share the 
same properties as our observations in terms of RV measurement 
uncertainties and time sampling, and that are drawn from known 
intrinsic distributions (see AppendixlBll. These tests also enabled 
us to look for systematics in our methods. From Table IBTI our 
approach tends to underestimate the intrinsic binary fraction by 
a couple of percentage points. The value obtained for the in¬ 
dex of the period distribution tt is typically 0.2 dex smaller than 
the adopted index. Similarly, the retrieved index of the mass ra¬ 
tio distribution k is typically 0.5 dex smaller than the adopted 
value. Nonetheless, these systematics remain smaller than the 
estimated statistical uncertainties on the measured value. 

We note that the k value estimated here differs to those fo und 
for Ga lactic binaries, e.g. k = — 0 .1 ± 0.6 (ISana et al.ll2012h and 
K'^0 (iKobulnickv et al.ll2014l). and k = —1.0 ± 0.4 for the O- 
type stars from the VFTS (IPaner Vllj) : this remains the case 
when taking into account the apparent bias to smaller values 
from our method (Appendix |B]l, as well as the omission of the 
five SB2 systems that will create an additional bias in the s ame 
sense. However, as discussed in Appendix C of IPaner Villi the 
Monte Carlo method is only weakly sensitive to the adopted 
value of K, with large uncertainties on its estimated value. 
Higher-quality spectroscopy (both in terms of S/N and cadence) 
of the B-type stars studied in this work would enable a more 
rigorous analysis of the mass-ratio distribution in the detected 
binaries. 

From these calculations, we estimate that the overall detec¬ 
tion probability of B-type binaries from the VFTS campaign is 
40 ± 10%. This suggests that we detected fewer than half the 
binaries with a B-type primary component. For the detected bi¬ 
nary fraction /B(obs) a detection probability p, a fraction 
[fobs X (1 -p)] / [p(l - fobs)] ^ 45% of these ‘single’ stars are 
expected to be binary. 


6. Summary 

We have investigated the multiplicity properties of 408 B-type 
stars in 30 Dor, employing the multi-epoch spectroscopy from 
the VFTS. We hnd that at least 25% of the stars display sig- 
niheant RV variations that are large enough to potentially be the 
spectroscopic signature of binarity. An additional 9% of the stars 
display RV variations - these are either binary candidates or dis¬ 
play intrinsic atmospheric variability. 

By modelling the properties of our observations, we estimate 
our observational biases and constrain the intrinsic multiplic¬ 
ity properties of the B-type stars in 30 Dor. The binary frac¬ 
tion of the unevolved (i.e. dwarf/giant) stars is estimated to be 
/B(true) =0.58 ±0.11, with a flat period distribution. The un¬ 
certainties on the mass-ratio distribution are large and the best- 
fit distribution should be considered with caution. Nonetheless, 
the multiplicity properties of the B-type stars agree (within the 
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Fig. 5. Projections of the global merit function. S', on the pairs of planes defined by tt, k and /B(true)- The location of the absolute 
maximum is indicated by a cross (+) and contours indicate loci of equal values in the merit function, expressed as a fraction of its 
absolute maximum (see inset in left-hand panel). 





0.1 1.0 10.0 100.0 
AHJD (day) 



0.1 1.0 10.0 100.0 
AHJD (doy) 


Fig. 6. Comparison of the observed (crosses) and simulated (lines) cumulative distributions of the peak-to-peak RV amplitudes 
(upper panels) and of the variability timescales (lower panels). The exponent tt (period distribution) is varied in the left-hand panels, 
while K (the exponent of the mass-ratio distribution), is varied in the right-hand panels. Values of tt, k and p are indicated in the 
upper-left insets of each panel. The lower-right insets indicate the VFTS detection probability for the range of parameters specified 
in Table |5] and for the adopted parent distributions. 














































































P. R. Dunstall et al.: Multiplicity of B-type stars in 30 Dor. 


Table 5. Properties of the Monte Carlo grid. The first three rows provide information on the orbital distributions: physical parameter 
(Col. 1), probability-density function (Col. 2) and applicability domain (Col. 3). Cols. 4 to 6 give the quantity, investigated range 
and step size in the grid. 


Parameter 

pdf 

Domain 

Var. 

Range 

Step 

P(d) 

(logioJ’)" 

0.15-3.5 

TT 

-2.50 - +2.50 

0.1 

g 

9" 

0.1 - 1.0 

K 

-3.50-+1.50 

0.1 

e 


10"® -0.9 

V 

-0.5 (fixed) 

nja 

/'B(true) 

nja 

n/a 

/'B(true) 

+0.40-+1.00 

0.02 


uncertainties) with those of the O-type stars in 30 Dor. Although 
the B-star sample is larger, the uncertainties are greater than for 
the O stars, a consequence of the B stars being fainter (hence 
having spectroscopy with generally lower S/N). Our results in¬ 
dicate that the important evolutionary implications of binarity 
inferred for the O-star sample probably also hold for the stars 
born as B-type, and hence for the entire population of supernova 
progenitors. Comprehensive spectroscopic monitoring of the B- 
type sample is now required to determine the period distribution 
to firmly test this conclusion. 
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Appendix A: Spatially-resolved companions 

The Medusa fibres subtend 1"2 on the sky which, at the distance 
of the LMC (50kpc), is equivalent to 0.29 pc. Some of the spec¬ 
tra could therefore be composites of both wide binaries and/or 
chance line-of-sight alignments. High-quality optical imaging 
from the wide-held F775W mosaic of 30 Dor taken with the 
Hubble Spac e Telescope (HST ) in programme GO-12499 (PI: 
Lennon; see ISabbi et 1201 3h includes 300 of the 403 stars 
(74%) with RV estimates. 

We checked the available shallow and deep images for com¬ 
panions within 1". Two of us (PLD and SdM) visually checked 
each of our targets for nearby companions which could have sig- 
nihcantly contaminated the LR02 spectroscopy; notable com¬ 
panions/features in this context are summarised in Table lA.ll 
Although subjective, the two independent evaluations were in 
excellent agreement, ensuring that important companions, as re¬ 
solved by HST imaging, are noted. 

We detect visual companions that could contribute to the ob¬ 
served spectra for thirty-hve of the B-star targets (i.e. 11.7 % of 
those with imaging available). In the context of searching for 
spectroscopic binaries such visual companions sample a very 
different range of p hysical separations (see e.g. Fig. 1 from 
ISana & Evansll20111) and they will not influence the results for 
the spectroscopic binary sample analysed in this paper. Even if 
they are physically bound, the RV amplitudes for these visual¬ 
binary systems will be below our detection threshold, and they 
will therefore be categorised as single using the spectroscopic 
criteria described in Section |3] (in the absence of other sources 
of RV variability). Such wide systems are unlikely to interact at 
any point in their evolution, and so will have no influence on 
our conclusions in this respect. However, to the extent that the 
observed spectra will be composite to some (generally small) 
degree, this contamination will need to be taken into account in 
any future atmospheric analyses of these spectra. 

Appendix B: Accuracy of the MC method 

As in iPaner Vlli we estimated the accuracy of the method by 
applying it to synthetic data which shared the same properties as 
the sample of early-B dwarfs analysed in this paper. We adopted 
different parent distributions to generate the synthetic observa¬ 
tions and investigated how well the method behaved in different 
parts of parameter space, including the region where the adopted 
parameters led to the best representation of the merit function. 
The results from these tests are reported in Table IB. Il 


Table A.l. B-type targets that may have a significant contribu¬ 
tion to their FLAMES spectra from nearby companions. 


VFTS 

Comment 

043 

Elongated image - close companion 

044 

Elongated image - close companion 

050 

Star with similar mag. at 

068 

Eainter star at ~0.5" 

127 

Eainter star at <0.5" 

133 

Star with similar mag. at ~0.5” 

167 

Elongated image - close companion 

212 

Elongated image - close companion 

238 

Two fainter stars, at <0.5''& ^0.5" 

276 

Star with similar mag. at ~0.5” 

278 

Crowded field, 4 fainter stars at <0.5” 

283 

Two stars with similar mag. at ~0.5” 

292 

Three fainter stars at <0.5” 

301 

Star with similar mag. at ~l” 

374 

Fainter star at <0.5” 

376 

Elongated image - close companion 

381 

Star with similar mag. at ~0.5” 

403 

Star with similar mag. at ~l” 

447 

Star with similar mag. at <0.5” 

448 

Star with similar mag. at <0.5” 

461 

Star with similar mag. at <0.5” 

463 

Fainter star at <0.5” 

480 

Two fainter stars at <0.5” 

504 

Star with similar mag. at ~l” 

548 

Fainter star at <0.5” 

553 

Fainter star at -^0.5” 

567 

Crowded field, star with similar mag. at 0.5” 

584 

Two fainter stars at ~0.5” 

632 

Fainter star at 0.5” 

643 

Star with similar mag. at ~0.5” 

644 

Star with similar mag. at ~0.5” 

662 

Fainter star at <0.5” 

712 

Star with similar mag. at <0.5” 

747 

Star with similar mag. at <0.5” 

833 

Elongated image - close companion 
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Table B.l. Overview of test results from synthetic datasets. The input multiplicity parameters are given first in Cols. 2 to 4. The 
medians and 0.16 and 0.74 percentiles of the retrieved parameters from a set of 50 test runs are then indicated. The properties of the 
computed Monte Carlo grid are the same as those in Table |5] 


ARVmin 

[kms-1] 

/B(true) 

Multiplicity properties 

TT 


16.0 

(/B(true) = 0.70) 
0.70 [ 0.62,0.82] 

(^Ttrue — 0.0) 

-0.20 [-0.50,+0.20] 

(^true — 0.0) 

-0.40 [-1.00,+0.90] 

16.0 

(/B(true) =0.70) 
0.68 [ 0.62,0.76] 

(^Ttrue — 0.5) 

-0.70 [-1.10,-0.40] 

(^true — 0.0) 

-0.50 [-1.00,+0.60] 

16.0 

(/B(true) =0.70) 
0.68 [ 0.54,0.78] 

(TTtrue = 0.0) 

-0.20 [-0.50,+0.20] 

(^true — 2.5) 

-2.90 [-3.50,-2.20] 

16.0 

(/B(true) =0.55) 
0.52 [ 0.42,0.64] 

(^Ttrue — 0.0) 

-0.20 [-0.60,+0.40] 

(^true — 2.5) 

-3.00 [-3.50,-2.00] 





